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Cell membranes interact via anchored receptor and ligand molecules. Central questions on cell adhesion 
concern the binding affinity of these membrane-anchored molecules, the mechanisms leading to the receptor- 
ligand domains observed during adhesion, and the role of cytoskeletal and other active processes. In this review, 
these questions are addressed from a theoretical perspective. We focus on models in which the membranes are 
described as elastic sheets, and the receptors and ligands as anchored molecules. In these models, the thermal 
membrane roughness on the nanometer scale leads to a cooperative binding of anchored receptor and ligand 
molecules, since the receptor-ligand binding smoothens out the membranes and facilitates the formation of 
additional bonds. Patterns of receptor domains observed in Monte Carlo simulations point towards a joint role 
of spontaneous and active processes in cell adhesion. The interactions mediated by the receptors and ligand 
molecules can be characterized by effective membrane adhesion potentials that depend on the concentrations 
and binding energies of the molecules. 



I. INTRODUCTION 

The adhesion of cells is mediated by the specific binding 
of receptor and ligand molecules anchored in the cell mem- 
branes. Cell adhesion processes are essential for the distinc- 
tion of self and foreign in immune responses, the formation of 
tissues, or the signal transduction across the synaptic cleft of 
neurons H]. These processes have therefore been sfiidied in- 
tensively with a variety of experimental methods ||2]|3]|4]|5]|6l. 
In addition, experiments on lipid vesicles with membrane- 
anchored receptor and ligand molecules aim to mimic the 
specific membrane binding processes leading to cell adhesion 

cma. 

In many adhesion processes, the anchored receptor and lig- 
and molecules can still diffuse, at least to some extent, within 
the contact area of the adhering membranes 13] HI |6l. As a 
consequence, the receptor-ligand complexes may form differ- 
ent spatial patterns such as clusters or extended domains in the 
contact area ll3ll6l [T0l[TTIl . These pattern formation processes 
can be understood in the framework of discrete models in 
which the membranes are divided into small patches, and the 
receptors and ligands are described as single molecules that 
are either present or absent in the patches 1 12l [T3lfT4l[T5l [T6^. 
These discrete models are lattice models on elastic surfaces, 
and have two advantages: (i) They automatically incorporate 
the mutual exclusion of receptor or ligand molecules anchored 
within the same membrane; and (ii) they lead to effective 
membrane adhesion potentials that provide an intuitive under- 
standing of the observed behavior in terms of nucleation and 
growth processes. 

Cell adhesion involves many different length scales, from 
nanometers to tens of micrometers. The largest length scales 
of micrometers correspond to the diameter of the cell and 
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the diameter of the contact zone in which the cell is bound 
to another cell or to a supported membrane. The separation 
of the two membranes in the cell contact zone is orders of 
magnitude smaller The membrane separation is comparable 
to the length of the receptor-ligand complexes, which have 
a typical extension between 15 and 40 nanometers |17|. Fi- 
nally, the smallest relevant length scale is the binding range 
of a receptor and a ligand molecule, i.e. the difference be- 
tween the smallest and the largest local membrane separation 
at which the molecules can bind. The binding range reflects 
(i) the range of the lock-and-key interaction, (ii) the flexibil- 
ity of the two binding partners, and (iii) the flexibility of the 
membrane anchoring. For the rather rigid protein receptors 
and ligands that typically mediate cell adhesion, the interac- 
tion range is around 1 nanometer. In contrast, the interaction 
range of surface-anchored flexible tether molecules with spe- 
cific binding sites is significantly larger 118., 19. .2011211 . 

The wide range of length scales has important conse- 
quences for modeling cell adhesion. In general, the elastic- 
ity of the cell membrane is affected by the bending rigidity k 
of the membrane 1.22,1 , the membrane tension cr, and the cy- 
toskeleton that is coupled to the membrane ll23l l24l l25l l26l . 
The tension dominates over the bending elasticity on length 
scales larger than the crossover length \J nja [27 1, which is of 
the order of several hundred nanometers for cell membranes 
ESl , while the bending elasticity dominates on length scales 
smaller than the crossover length. The elastic contribution of 
the cytoskeleton is relevant on length scales larger than the av- 
erage distance of the cytoskeletal membrane anchors, which 
is around 100 nanometers The overall shape of the cell 
membrane on micrometer scales therefore is governed by the 
cytoskeletal elasticity and membrane tension. In the cell ad- 
hesion zone, however, the relevant shape deformations and 
fluctuations of the membranes occur on length scales up to 
the average distance of the receptor-ligand bonds, since the 
bonds locally constrain the membrane separation. The aver- 
age distance of the bonds roughly varies between 50 and 100 
nanometers for typical bond concentrations in cell adhesion 
zones |3|. The relevant membrane shape deformations and 
fluctuations in the cell contact zone are therefore dominated 
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by the bending rigidity of the membranes. 

The adhesion of cells is mediated by a multitude of differ- 
ent receptor and ligand molecules. Some of these molecules 
can be strongly coupled to the cytoskeleton. In focal adhe- 
sions of cells, for example, clusters of integrin molecules are 
tightly coupled to the cytoskeleton via supramolecular assem- 
blies that impose constraints on the lateral separation of the 
integrins [5 29 1. Through focal adhesions, cells exert and 
sense forces IMl EH EH ifS] |3l ED . Other receptor and lig- 
and molecules are not or only weakly f36l coupled to the 
cytoskeleton. These molecules are mobile and diffuse within 
the membranes. The diffusion process can be observed with 
single-molecule experiments |37u38J . 

The adhesion of membranes via mobile receptor and lig- 
and molecules has been studied theoretically with a variety 
of models. These models can be grouped into two classes. 
In both classes, the membranes are described as thin elas- 
tic sheets. In the first class of models, the description is 
continuous in space, and the distribution of the membrane- 
anchored receptor and ligand molecules on the membranes 
are given by continuous concentration profiles Il39ll40l 14111421 
[43 44 45 46 1. Dynamic, time-dependent properties of such 
models have been studied by numerical solution of reaction- 
diffusion equations 1451 l47l l48l |49J . In the second, more 
recent class of models, the membranes are discretized, and 
the receptors and ligands are described as single molecules 
lIiainillSlIISIiaEOllSllEailllEllia The dynamic 
properties can be numerically studied with Monte Carlo sim- 
ulations lfT4ll28ll53ll54l [56l. and central aspects of the equi- 
librium behavior can be directly inferred from the partition 
function lfT3l [Bl [161 l50l . 



II. EFFECTIVE ADHESION POTENTIAL 
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FIG. 1: (a) A membrane segment with receptor molecules (top) 
interacting with ligands embedded in an apposing membrane (bot- 
tom). A receptor can bind a ligand molecule if the local separation 
of the membranes is close to the length of the receptor-ligand com- 
plex. - (b) The attractive interactions between the receptor and lig- 
and molecules lead to an effective single-well adhesion potential V^i 
of the membranes. The depth Ud of the potential well depends on 
the concentrations and binding affinity of receptors and ligands, see 
eq. |T]). The width /„c of the binding well is equal to the binding 
range of the receptor-ligand interaction. 



In discrete models, the two apposing membranes in the con- 
tact zone of cells or vesicles are divided into small patches 
|[I2[l3[l5][IS|28l|50l|5l]|52l|53l|53|5^ These patches can 
contain a single receptor or ligand molecule. Mobile recep- 
tor and ligand molecules diffuse by 'hopping' from patch to 
patch, and the thermal fluctuations of the membranes are re- 
flected in variations of the local separation of apposing mem- 
brane patches. A receptor can bind to a ligand molecule if the 
ligand is located in the membrane patch apposing the recep- 
tor, and if the local separation of the membranes is close to 
the length of the receptor-ligand complex. In these models, 
the linear size a of the membrane patches is typically chosen 
around 5 nm to capture the whole spectrum of bending defor- 
mations of the lipid membranes |57|. 

Cells can interact via a multitude of different receptors and 
ligands. However, it is instructive to start with the relatively 
simple situation in which the adhesion is mediated by a single 
type of receptor-ligand bonds as in fig. [TJ^)- Such a situation 
occurs if a cell adheres to a supported membane with a single 
type of ligands, or if a vesicle with membrane-anchored re- 
ceptors adheres to a membrane with complementary ligands. 
The effective membrane adhesion potential mediated by the 
receptor and ligands can be calculated by integrating over all 



possible positions of the receptors and ligands in the partition 
function of the models ||T6ll50l . In the case of a single type 
of receptors and ligands, the effective adhesion potential of 
the membranes is a single-well potential with the same range 
Zwe as the receptor-ligand interaction, but an effective binding 
energy [/gf that depends on the concentrations and binding en- 
ergy U of receptors and ligands, see fig.[TJb). For typical con- 
centrations of receptors and ligands in cell membranes, which 
are more than two orders of magnitude smaller than the max- 
imum concentration ~ 4 • 10^//im^ in our discretized 
membranes with patch size a ~ 5 nm, the effective potential 
depth is 1 58 1 

where [i?] and [L] are the area concentrations of unbound re- 
ceptors and ligands. The quantity 

ifpi = a^e^/^-'^^ (2) 

in eq. ([T]) can be interpreted as the binding equilibrium con- 
stant of the receptors and ligands in the case of two planar and 
parallel membranes with a separation equal to the length of the 
receptor-ligand bonds. The equilibrium constant characterizes 
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FIG. 2: (a) Two membranes interacting via long (red) and short 
(green) receptor-ligand complexes. - (b) The attractive interactions 
between the two types of receptors and ligands lead to an effective 
double-well adhesion potential Kf of the membranes. The poten- 
tial well 1 at small membrane separations I reflects the interactions 
of the short receptor-ligand complexes, and the potential well 2 at 
larger membrane separations the interactions of the long receptor- 
ligand complexes. The depths Uf and U2 of the two potential wells 
depend on the concentrations and binding energies of the two types 
of receptors and ligands, see eqs. ^ and |4|. 



the binding affinity of the molecules and can, in principle, be 
measured with the surface force apparatus in which the ap- 
posing membranes are supported on rigid substrates |59, 60|. 
In the case of flexible membranes, the binding affinity of the 
receptors and ligands is more difficult to capture, see next sec- 
tion. 

The interaction of cells is often mediated by several types 
of receptor-ligand complexes that differ in their length. For 
two types of receptors and ligands as in fig. |2] the effective 
adhesion potential of the membranes is a double-well poten- 
tial lil5J . The depths of the two wells 



(3) 
(4) 



depend on the concentrations and binding energies Ui and 
U2 of the different types of receptors and ligands |58|. In 
analogy to eq. Ol, the quantities Kp\ i = a^^^i/^BT 
Kp\ 2 = o?gU2/kBT ^^jj interpreted as binding equilibrium 
constants in the case of planar membranes with a separation 
equal to the lengths li or I2 of the receptor-ligand complexes. 

Repulsive membrane-anchored molecules such as anchored 
polymers or glycoproteins can lead to additional barriers in 
the effective adhesion potential ifTSl [50l ISTl . The effective 
adhesion potentials simplify the characterization of the equi- 



librium properties of the membranes, and lead to an intuitive 
understanding of these properties, see next sections. 



ni. BINDING COOPERATIVITY 

A receptor molecule can only bind an apposing ligand if 
the local membrane separation is comparable to the length of 
the receptor-ligand complex. A central quantity therefore is 
the fraction P\j of the apposing membranes with a separation 
within the binding range of the receptor-ligand interaction. 
The concentration of bound receptor-ligand complexes 



[RL] 



PbKpi [R][L] 



(5) 



is proportional to as well as to the concentrations [R] and 
[L] of unbound receptors and ligands |58 1. 

Thermal shape fluctuations of the membranes on nanome- 
ter scales in general lead to values of Pi, smaller than 1 . For 
cell membranes, these nanometer scale fluctuations are not, 
or only weakly, suppressed by the cell cytoskeleton, in con- 
trast to large-scale shape fluctuations. For simplicity, we as- 
sume now that the adhesion of the membranes is mediated by 
a single type of receptor and ligand molecules as in fig. [TJa)- 
The precise value of Pb then depends on the well depth Ugf 
of the effective adhesion potential shown in fig. [ijb), and on 
the bending rigidities of the membranes. For typical lengths 
and concentrations of receptors and ligands in cell adhesion 
zones, the fraction Pf, of the membranes within binding range 
of the receptors and ligands turns out to be much smaller than 
1, and scaling analysis and Monte Carlo simulations lead to 
the relation 



Pb-ci^ll^U,i/{kBT)^ (6) 
with prefactor c = 13 ± 1. With eqs. ([T]i and (j2|, we obtain 



Pb^c{n/kBT)ll^Kp, [R][L] 



(7) 



which shows that the membrane fraction P^ within the binding 
range of the receptors and ligands is proportional to \R\ and 
[L] . Inserting eq. (j7| into eq. (jsjl leads to 



[RL]^c{K/kBT)ll^Kl[Rf[L] 



(8) 



The concentration [RL] of receptor-ligand complexes in the 
adhesion zone thus depends quadratically on the concentra- 
tions [R\ and [L] of unbound receptors and ligands, which 
indicates cooperative binding. The binding cooperativity re- 
sults from a 'smoothening' of the thermally rough membranes 
and, thus, an increase of Pi, with increasing concentrations \R\ 
and [L] of receptors and ligands, which facilitates the forma- 
tion of additional receptor-ligand complexes. The relations 
(|7ji and (jsj are good approximations up to Pb ^ 0.2, and can 
be extended to larger values of Pb ll58l . 

For soluble receptor and ligand molecules, in contrast, the 
volume concentration of the bound receptor-ligand complexes 



[i?L]3D = i^3D[i?]3D[i]3D 



(9) 
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is proportional to the volume concentrations [R]iD and [L]jo 
of unbound receptors and unbound ligands in the solution. 
The binding affinity of the molecules then can be character- 
ized by the equilibrium constant K^r,, which depends on the 
binding free energy of the complex 1,61, 62 63 1. In analogy to 
eq. the binding affinity of membrane-anchored receptors 
and ligands is often quantified by 



K2D = 



[RL] 



mi] 



(10) 



where [RL], [R], and [L] are the area concentrations of bound 
receptor-ligand complexes, unbound receptors, and unbound 
ligands 1641 l65l l66l . 

However, it follows from our relation ([SJ that K2D is not 
constant, but depends on the concentrations of the receptors 
and ligands. From the eqs. Q and (j5]l, we obtain the general 
relation 



K 



2D 



PhK. 



(11) 



As mentioned in the previous section, Kp\ is the well-defined 
two-dimensional equilibrium constant of the receptors and lig- 
ands in the case of planar membranes with = 1. e.g. two 
supported membranes in the surface force apparatus with a 
separation equal to the length of the receptor-ligand complex 



The relation ( [TT| ) also helps to understand why different ex- 
perimental methods for measuring K2D have led to values that 
differ by several orders of magnitude ll65l . In fluorescence re- 
covery experiments, K2D is measured in the equilibrated con- 
tact zone of a cell adhering to a supported membrane with 
fluorescently labeled ligands 1671 l68l l69l iTOl . In micropipette 
experiments, in contrast, K2D is measured for initial contacts 
between two cells |66, 71, 72 1, for which Pf, can be several 
orders of magnitude smaller than in equilibrium liSSl . 



IV. ADHESION OF VESICLES 

Important aspects of cell adhesion can be mimicked by Upid 
vesicles with anchored receptor molecules ||7] [8] |9j |73] |74| 
l75l l76l . We focus here on a vesicle adhering to a supported 
membrane with complementary ligands. In the strong adhe- 
sion limit, the shape of the vesicle can be approximated by 
a spherical cap 1771 ITSl . The volume of the cap depends on 
the osmotic pressure balance between the outside and the in- 
terior of the vesicle. If this volume is nearly constant, the con- 
tact area Ac is nearly independent of the adhesion free energy 
Il52l|78 79 80|. 

Since the total number Nb. of receptors in the vesicle mem- 
brane is fixed, we have 



Nr^[R]A+[RL]Ac 



(12) 



where A is the total area of the vesicle, and [RL] is the concen- 
tration of receptor-ligand complexes in the contact area. For 
typical small concentrations of receptors and ligands, the con- 
centration [R] of unbound receptors within the contact area 
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FIG. 3: Concentrations [RL] and [7?] of bound and unbound recep- 
tors as a function of the receptor number Nr of a vesicle adhering 
to a supported membrane, see eqs. |T3| and ( |14^ . In this numerical 
example, we have chosen the bending rigidity k — 25fcflr, the bind- 
ing range l„c = 1 nm, the binding equilibrium constant Kpi = 1 fim^ 
for planar membranes, and a ligand concentration [L] — 20//im^ for 
the supported membrane, which results in the value 0.14//im^ for 
the parameter b in eqs. \13\ and \14\ . The total area of the vesicle is 
100 fim^, and the area of the contact zone is 30 /im^. The fraction Pi, 
of the vesicle membrane in the contact zone with a separation within 
binding range of the receptors and ligands varies with [R], see eq.jT] 
and attains the maximum value Pb = 0.17 for Nr — 5000 in this 
example. 



and within the non-adhering membrane section of the vesi- 
cle are approximately equal in equilibrium since the excluded 
volume of the receptor-ligand complexes in the contact area is 
negligible. With eq. ([8]l, we then obtain 



and 



[R] 



[RL] 



^JA^ + 4Mc7Vfl - A 
26X 



(VA2 + 46A,7Vfl,- A)^ 
46l2 



(13) 



(14) 



with 6 = c{K.lkBT)l\^K'^^V^ . Here, [U] is the concentration 
of the unbound ligands in the supported membrane, which is 
nearly independent of the binding state of the vesicle if the 
membrane is large. Because of the binding cooperativity of 
the receptors and ligands, the concentrations [i?] and [RL] are 
not hnear in iV/{, see fig.[3]for a numerical example. 



V. DOMAINS OF LONG AND SHORT 
RECEPTOR-LIGAND COMPLEXES 

A. Critical concentrations for domain formation 

Cells often interact via receptor-ligand complexes that dif- 
fer significantly in size. For example, two important com- 
plexes in T-cell adhesion are the complexes of the T-cell re- 
ceptor (TCR) with a length of 15 nm and integrin complexes 
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FIG. 4: 



Phase diagram of membranes adhering via long and short receptor/hgand complexes. 



The membranes are unbound for small well 



3cept( 

depths Ui and f/? of the effective interaction potential shown in Fig.bfb), i.e. for small concentrations or binding energies of receptors and 
ligands, see eqs. ([sjl and |4|. At large values of Ui and (7|'> the membranes are either bound in well 1 or well 2, i.e. they are either bound by 
the short or by the long receptor/ligand complexes. At intermediate well depths Uf and U2, the membranes are bound in both potential wells. 
The critical point for the lateral phase separation (star) follows from eq. ^15) . For typical dimensions of cell receptors and ligands, the critical 
well depth for lateral phase separation is significantly larger than the critical depths of unbinding 1151 . In the absence of other repulsive 
interactions as assumed here, the membranes unbind due to steric repulsion. 



with a length of 40 nm ifTTl . The length mismatch induces 
a membrane-mediated repulsion between the different com- 
plexes because the membranes have to be curved to compen- 
sate the mismatch, which costs bending energy. 

The equilibrium behavior of two membranes adhering via 
long and short receptor-ligand complexes is determined by the 
effective double-well adhesion potential shown in fig. |2|b), 
and by the bending rigidities of the membranes. The depths 
of the two wells reflect the concentrations and binding affin- 
ity of the two different types of receptors and ligands, see 
eqs. ([3]) and Q in section |ll] If the two wells are relatively 
shallow, membrane segments can easily cross the barrier be- 
tween the wells, driven by thermal fluctuations. If the two 
wells are deep, the crossing from one well to the other well 
is suppressed by the potential barrier of width 4a between the 
wells. The potential barrier induces a line tension between a 
membrane segment that is bound in one of the wells and an 
adjacent membrane segment bound in the other well. Beyond 
a critical depth of the potential wells, the line tension leads 
to the formation of large membrane domains that are bound 
in well one or well two, see fig. [4] Within each domain, the 
adhesion of the membranes is predominantly mediated by one 
of the two types of receptor-ligand complexes. 

Scaling arguments indicate that domains bound in either 
well 1 or well 2 are formed if the depths Uf and C/|' of the 
two potential wells exceed the critical potential depth ifTSl 



cjkBT)^ 



(15) 



Numerical results from Monte Carlo simulations confirm 
eq. (fTSl) and lead to the value c = 0.225 ± 0.02 for the di- 



mensionless prefactor ifTSl . The critical potential depth thus 



depends on the effective rigidity k = kiK2/{ki + K2) of 
two membranes with bending rigidities ki and K2 and the 
width and separation 4a of the two potential wells, see 
fig.|2|b). The separation /ba of the wells is close to the length 
mismatch of the different types of receptor-ligand complexes, 
which is 25 nm for T cells. A reasonable estimate for the in- 
teraction range Z„e of the protein receptors and ligands that 
mediate cell adhesion is 1 nm, see section |l] With an effec- 
tive bending rigidity k = kiK2/{ki + K2) of, e.g., 25 /csT, 
we obtain the estimate Uf ~ 360 kBT/ firn^ for the criti- 
cal potential depth of domain formation during T-cell adhe- 
sion. For planar-membrane equilibrium constants K^u and 
Kp\ 2 around 1 firn^, for example, the effective potential depths 
Q and (|4]) exceed this critical potential depth if the con- 
centrations of unbound receptors and ligands are larger than 
20//im2. 



B. Domain patterns during immune cell adhesion 

The domains of long and short receptor-ligand complexes 
formed during the adhesion of T cells and other immune cells 
such as natural killer cells evolve in characteristic patterns. 
For T cells, the domains either contain complexes of TCR and 
its ligand MHC-peptide, or integrin complexes. The final do- 
main pattern in the T-cell contact zone is formed within 15 to 
30 minutes and consists of a central TCR domain surrounded 
by a ring-shaped integrin domain |3, 10 1, see fig.|5ja). Inter- 
estingly, the intermediate patterns formed within the first min- 
utes of T-cell adhesion are quite different 13] [lT|. They are 
either inverse to the final pattern, with a central integrin do- 
main surrounded by a ring-shaped TCR domain, see fig. |5jc), 
or exhibit several nearly circular TCR domains in the contact 
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hesion, the Monte Carlo patterns exhibit small TCR clusters 
llT4ll . In the following seconds, the diffusion of free TCR and 
MHC-peptide molecules into the contact zone lead to faster 
growth of TCR clusters close to the periphery of the contact 
zone [84|. For sufficiently large TCR-MHC-peptide concen- 
trations, the peripheral TCR clusters grow into the ring-shaped 
domain of fig. |5|c). At smaller TCR-MHC-peptide concen- 
tration, the initial clusters evolve into the pattern of fig.|5jd). 
In agreement with experimental observations 111], only these 
two types of intermediate patterns are formed in the simula- 
tions. The simulated patterns emerge spontaneously from the 
nucleation of TCR clusters and the diffusion of unbound TCR 
and MHC-peptide into the contact zone. 



FIG. 5: Domain patterns in the T-cell contact zone: (a) Final pattern 
of helper T cells with a central TCR domain (green) surrounded by an 
integrin domain (red) l3l 1101 . The pattern results from cytoskeletal 
transport of TCRs towards the contact zone center l6lll4l . - (b) Sim- 
ulated final pattern in the absence of TCR transport 1141 . The length 
of the boundary line between the TCR and the integrin domain is 
minimal in this pattern. - (c) and (d) The two types of intermediate 
patterns observed in the first minutes of adhesion II II . In simula- 
tions, both patterns result from the nucleation of TCR clusters in the 
first seconds of adhesion and the subsequent diffusion of unbound 
TCR and MHC-peptide ligands in the contact zone 1 14|. The closed 
TCR ring in pattern (c) forms from fast-growing TCR clusters in the 
periphery of the contact zone at sufficiently large TCR-MHC-peptide 
concentrations. The pattern (d) forms at smaller TCR-MHC-peptide 
concentrations. 



zone, see fig.[5]^d). 

To understand these patterns, several groups have modeled 
and simulated the adhesion of T cells and other immune cells 
ni SI ill ESI |53] |54l |56l 11] 111. One open question 
concerned the role of the T-cell cytoskeleton, which polarizes 
during adhesion, with a focal point in the center of the con- 
tact zone I r,'83l. Some groups have found that the final T-cell 
pattern with a central TCR domain can emerge independently 
of cytoskeletal processes BtI IsTI . In contrast, Monte Carlo 
simulations of discrete models indicate that the central TCR 
cluster is only formed if TCR molecules are actively trans- 
ported by the cytoskeleton towards the center of the contact 
zone L14J . The active transport has been simulated by a bi- 
ased diffusion of TCRs towards the contact zone center, which 
implies a weak coupling of TCRs to the cytoskeleton. In the 
absence of active TCR transport, the Monte Carlo simulations 
lead to the final, equilibrium pattern shown in fig. |5|b), which 
minimizes the energy of the boundary line between the TCR 
and the integrin domain 1 14|. In agreement with these simu- 
lations, recent T-cell adhesion experiments on patterned sub- 
strates reveal cytoskeletal forces that drive the TCRs towards 
the center of the contact zone |l6l[36l. The experiments indi- 
cate a weak frictional coupling of the TCRs to the cytoskeletal 
flow I.36J. 



The intermediate patterns formed in the Monte Carlo sim- 
ulations closely resemble the intermediate immune-cell pat- 
terns shown in figs. |5|c) and (d). In the first seconds of ad- 



C. Implications for T-cell activation 

T cells mediate immune responses by adhering to cells that 
display foreign peptides on their surfaces HI |85|. The pep- 
tides are presented by MHC proteins on the cell surfaces, and 
are recognized by the T-cell receptors (TCR). T-cell activa- 
tion requires the binding of TCRs to the MHC-peptide com- 
plexes. But how precisely these binding events trigger T-cell 
activation still is a current focus of immunology (for reviews, 
see refs. [861 l87l [88l ). Recent experiments indicate that the 
first T-cell activation signals coincide with the formation of 
TCR microclusters within the first seconds of T-cell adhesion 
189,90 91 92,93,94|. 

In the discrete model introduced in section [11] and fig. [2] 
TCR clusters in the T-cell contact zone can only form if two 
conditions are met. First, the effective potential depth Uf for 
the short TCR-MHC-peptide complexes and the depth C/j* for 
the long integrin complexes have to exceed the critical depth 



formation of TCR clusters 
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FIG. 6: Schematic diagram for the joint role of foreign and self 
MHC-peptides in TCR cluster formation and T-cell activation. Here, 
[Alf] is the concentration of foreign MHC-peptides, and [AIs] is the 
concentration of self MHC-peptides. The solid line represents the 
threshold for TCR cluster formation given by eq. (jT7|. The slope 
of this line is the negative ratio Kjus/Kjuf of the binding equilib- 
rium constants for the interaction of TCR with self MHC-peptide and 
with foreign MHC-peptide. For simplicity, we have assumed here a 
single, dominant type of self MHC-peptides. 
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(15 I. Second, the effective potential Uf for the TCR com- 
plexes has to be larger than the effective depth U2 in the sit- 
uation where no TCRs are bound. To understand the second 
condition, one has to realize that the concentrations of un- 
bound TCRs and unbound integrins depend on the area frac- 
tions of the TCR and integrin clusters and domains in the con- 
tact zone. If no TCRs are bound, i.e. if the whole contact zone 
is occupied by an integrin domain, the concentration of un- 
bound TCRs is maximal. Hence, also the effective depth Uf 
for the TCRs is maximal in this situation, see eq. Q. TCR 
clusters now form if Uf is larger than C/|^, which leads to a 
decrease in the concentration of unbound TCRs and, thus, to a 
decrease in Uf. The area fraction of the TCR clusters grows 
until the equilibrium situation with Uf = Uf is reached ifTSl . 

T-cell activation requires a threshold concentration of for- 
eign MHC -peptide complexes. Interestingly, the threshold 
concentration of foreign MHC-peptide depends on the con- 
centration of self MHC-peptide complexes, i.e. of complexes 
between MHC and self peptides derived from proteins of the 
host cell II95I |96]| . The foreign MHC-peptide complexes, in 
contrast, are complexes of MHC with peptides derived from 
viral or bacterial proteins. Self MHC-peptides typically bind 
weakly to TCR, since strong binding can result in autoimmune 
reactions. However, the number of self MHC-peptide com- 
plexes typically greatly exceed the number of foreign MHC- 
peptide on cell surfaces. Both self and foreign MHC-peptide 
complexes contribute to the effective potential depth Uf of 
the TCR-MHC-peptide interaction. For simplicity, we assume 
here a single, dominant type of self MHC-peptides with con- 
centrations [Mg\. The effective potential depth then is 

Uf = kBT{[T][Mf]K^Mf + [T][M,]KjMs + ■ ■ ■) (16) 

where [T] and [Mf] are the concentrations of unbound TCR 
and foreign MHC-peptide, and A'xMf and /^tms are the bind- 
ing equilibrium constants of foreign and self TCR-MHC- 
peptide complexes in the case of planar membranes, see 
The dots in eq. 



16 1 indicate possible contributions 



eq 

from other receptor-ligand complexes with the same length 
as the TCR-MHC-peptide complex, e.g. from the CD2-CD58 
complex ifTTl . In addition, repulsive glycoproteins with a 
length larger than TCR-MHC-peptide complex can affect Uf 
lfT4ll56l . Similarly, the depth Uf of the second well depends 
on the concentrations and binding equiUbrium constants of in- 
tegrins and its ligands. 

Let us now suppose that the numbers of TCRs, co-receptors 
such as CD2, integrins, and glycoproteins are approximately 
equal for different T cells and apposing cells, while the num- 
bers of foreign and self MHC-peptides vary. The second con- 
dition uf > uf for TCR cluster formation then leads to 



(17) 



where C( is a dimensionless threshold that depends on the 
TCR, co-receptor, integrin, and glycoprotein concentrations, 
etc. The threshold concentration of foreign MHC-peptide 
complexes for TCR cluster formation thus depends on the 
concentration of self MHC-peptide complexes, see fig. [6] If 
the formation of TCR microclusters coincides with early ac- 
tivation signals as suggested in Refs. |89l |90ll , the inequality 



(17 1 also helps to understand the joint role of foreign and self 
MHC-peptides in T-cell activation. 



VI. ADHESION VIA CROSSLINKER MOLECULES OR 
ADSORBED PARTICLES 

The binding of receptor molecules on apposing membranes 
or surfaces is sometimes mediated by linker or connector 
molecules, see fig.|7|a). Biotinylated lipids in apposing mem- 
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FIG. 7: (a) Two membranes with receptors binding to solute 
molecules. At small membrane separations, the molecules can bind 
to two apposing receptors and, thus, crosslink the membranes. - (b) 
Two membranes in contact with a solution of adhesive molecules or 
particles. A particle can bind the two membranes together for mem- 
brane separations slightly larger than the particle diameter. At larger 
separations, the particles can only bind to one of the membranes. 
- (c) Effective adhesion potential V of the membranes in subfigure 
(b) as a function of the membrane separation I for small concentra- 
tions of the particles |101| . The effective potential has a minimum 
at the separation I = d + r where d is the particle diameter, and 
r is the range of the adhesive interaction between the particles and 
the membranes. At this separation, the particles are bound to both 
membranes. The effective potential is constant for large separations 
at which the particles can only bind to one of the membranes. The 
potential barrier at intermediate separations d-\-2r <l <2d results 
from the fact that a particle bound to one of the membranes locally 
'blocks' the binding of other particles to the apposing membrane. 
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FIG. 8: Effective adhesion energy Ud, given in eq. \20\ , as a function of the volume fraction of the adhesive particles for the binding energy 
U — SksT and q = 0.25. The effective adhesion energy is maximal at the optimal volume fraction (p* ~ ^^^"^ j q ~ 1.34 • 10~^. At the 
optimal volume fraction, the particle coverage of two planar parallel membranes is close to 50% for large separations, and almost 100% for 
small separations at which the particles can bind to both surfaces llOll . 



branes, for example, can be crosslinked by the connector 
molecule streptavidin ||2l |92l . The effective binding affinity 
of the membranes then depends both on the area concentra- 
tions of the membrane receptors and the volume concentra- 
tion of the linker molecules. A similar situation arises if ad- 
hesive molecules or particles directly bind to lipid bilayers 
(98 99, 100|. The adhesive particles can crosslink two appos- 
ing membranes if the membrane separation is close to the par- 
ticle diameter, see fig. |7|b). At large membrane separations, 
the particles can only bind to one of the membranes. 

The effective, particle-mediated adhesion potential of the 
membranes can be determined by integrating over all possi- 
ble positions of the adhesive particles or linker molecules in 
the partition function of the considered model. Conceptually, 
this is similar to the calculation of the effective adhesion po- 
tential for membranes interacting via anchored receptors and 
ligands, which requires an integration over all positions of the 
receptor and ligand molecules in the membranes, see section 
[n] For simplicity, we consider here the adhesive particles of 
fig. |7]^b), which interact directly with the lipid bilayers. The 
explicit integration over the particle positions requires spatial 
discretizations. In a lattice model, the space between the ap- 
posing membranes is discretized into a cubic lattice with a lat- 
tice spacing equal to the particle diameter d IIIOII . In an alter- 
native semi-continuous model, only the two spatial directions 
parallel to the membranes are discretized, while the third spa- 
tial direction perpendicular to the membranes is continuous 
mOlll . In both models, the effective, particle-mediated adhe- 
sion potential at large membrane separations has the form 



1 



,U/kBT 



(18) 



for small volume fractions (j) and large binding energies U of 
the particles. At small separations close to particle diameter, 
the adhesion potential exhibits a minimum 



dp 



■In 



1 



(19) 



The model-dependent factor q in eqs. (18 1 and ( [T9| ) has the 
value 1 in the lattice gas model and the value r/d in the semi- 
continous model with interaction range r of the adhesive par- 
ticles. In the semi-continuous model, the potential minimum 
is located at the membrane separation ^ = d + r, see figjTlc). 

Independent of these two models, the eqs. ( [T8] l and ([T^jican 
also be understood as Langmuir adsorption free energies per 
binding site. Eq. ( 19 1 can be interpreted as the Langmuir ad- 



sorption free energy for small membrane separations at which 
a particle binds both membranes with total binding energy 
1U , and eq. ( 18 1 as the Langmuir adsorption free energy for 



large separations. The Langmuir adsorption free energies re- 
sult from a simple two-state model in which a particle is either 
absent (Boltzmann weight 1) or present (Boltzmann weights 
q (f> e^v/ksT ^jjj ^ ^ ^u/ksT ^ respectively) at a given binding 
site. 

The effective, particle-mediated adhesion energy of the 
membranes can be defined as 



In 



{l + qcj)eU/kBTy 



(20) 



Interestingly, the effective adhesion energy is maximal at 

/q, and considerably 



the volume fraction 



-U/ksT 



smaller at smaller or larger volume fractions, see fig. [8] At 
this optimal volume fraction, the particle coverage Coc — 
-{(P/2){dV^/dU) « + of the unbound mem- 

branes is 50%. In contrast, the particle coverage Cmin = 
-{d^/2){dVn,iJdU) w </./(0 + (j)*e-^/''^'^) of the bound 
membranes is close to 100% at = 0*. Bringing the surfaces 
from large separations within binding separations thus does 
not 'require' desorption or adsorption of particles at the opti- 
mal volume fraction. The existence of an optimal particle vol- 
ume fraction has important implications that are accessible to 
experiments, such as 're-entrant transitions' in which surfaces 
or colloidal objects first bind with increasing concentration of 
adhesive particles, and unbind again when the concentration 
is further increased beyond the optimal concentration. 



VII. ACTIVE SWITCHING OF ADHESION RECEPTORS 



Some adhesion receptors can be switched between dif- 
ferent conformations. A biological example of switchable, 
membrane-anchored adhesion receptors are integrins. In one 
of their conformations, the integrin molecules are extended 
and can bind to apposing ligands [i l02lll03lll04L In another 
conformation, the molecules are bent and, thus, deactivated. 
The transitions between these conformations are triggered by 
signaling cascades in the cells, which typically require energy 
input, e.g. via ATP. Because of this energy input, the switching 
process is an active, non-equilibrium process. In biomimetic 
applications with designed molecules, active conformational 
transitions may also be triggered by light |i,105..106J . Other ac- 
tive processes of biomembranes include the forces exerted by 
embedded ion pumps |[l02l[I08l|IMl[Iini[rrT'Tl2l or by the 
cell cytoskeleton | iTJl ll 141 . see also section V B Active con- 
formational transitions of membrane proteins have also been 
suggested to couple to the local thickness 111 151 or curvature 
fl 16,1 of the membranes. 



In the absence of active processes, the adhesiveness of 
two membranes with complementary receptor and ligand 
molecules depends on the concentration and binding ener- 
gies of the molecules, and can be captured by effective adhe- 
sion potentials, see section|ll] The adhesiveness of membranes 
with actively switched receptors, in contrast, depends also on 
the switching rates of the receptors. In the example illustrated 
in fig.[9j the anchored receptors can be switched between two 
states: an extended 'on' -state in which the receptors can bind 
to ligands anchored in the apposing membrane, and an 'off- 
state in which the receptors can't bind. In this example, the 
switching process from the on- to the off-state requires en- 
ergy input, e.g. in the form of ATP. As a consequence, the rate 
for this process does not depend on wether the receptor 
is bound or not, in contrast to equilibrium situations without 
energy input. In an equilibrium situation, the rate for the tran- 
sition from the on- to the off-state depends on the binding state 
and binding energy of a receptor 



The active switching of the receptors enhances the shape 
fluctuations of the membranes 111 171 11181 11 191 . Since the 
steric repulsion of the membranes increases with the shape 
fluctuations, this enhancement of shape fluctuations leads to 
larger membrane separations. In Fig.|9jb), the average mem- 
brane separation is shown as a function of the switching rates 
for equal on- and off -rates aj+ = In this example, the 
fractions of receptors in the on- and off-state are constant and 
equal to 0.5. The active switching leads to a stochastic reso- 
nance of the membrane shape fluctuations, with a maximum 
of the membrane separation at intermediate switching rates. 
At the resonance point, the switching rates are close to the 
fluctuation relaxation rate l/r of a membrane segment with 
a linear size equal to the average separation of the receptors 
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FIG. 9: (Top) A membrane with switchable receptors adhering to 
a second membrane with complementary ligands. The receptors 
are switched between a stretched, active conformation and a bent, 
inactive conformation. In the stretched conformation, the adhe- 
sion molecules can bind to their ligands in the apposing membrane. 
" (Bottom) Monte Carlo data for the average reseated membrane 
separation z = T/a^J K/knT as a function of the switching rate 
uj = ujj^ — tj_ of the receptors. Here, lj+ and tj_ are the on- 
and off-switching rates of the receptors, and r is the characteristic 
relaxation time of a membrane segment with a linear size equal to 
the mean distance of the receptors. The active switching leads to a 
stochastic resonance with increased membrane separations at inter- 
mediate switching rates. The details of the Monte Carlo simulations 
are described in Ref. II 171 . In this example, the binding energy of 
the receptors and ligands isU — 2.8kBT. 



VIII. DISCUSSION AND CONCLUSIONS 

We have reviewed theoretical models for the adhesion of 
biomimetic membranes and cells via anchored but mobile re- 
ceptor and ligand molecules. In these models, the membranes 
are described as elastic surfaces, and the receptors and ligands 
as single molecules. We have argued in the introduction that 
the elasticity of the membranes is dominated by their bend- 
ing energy on the relevant lateral length scales up to average 
separation of the receptor-ligand complexes, which is between 
50 and 100 nm for typical concentrations of the complexes in 
cell adhesion zones [3J. The crossover length y^^/'''' above 
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which the tension a dominates over the bending rigidity k, 
is clearly larger for typical membrane tensions cr, see intro- 
duction. However, the average separation of cytoskeletal an- 
chors in cell membranes may be close to the average separa- 
tion of the receptor-ligand complexes. In the absence of active 
processes, the coupling of the membrane to the cytoskeleton 
may lead to a suppression of membrane shape fluctuations on 
length scales larger than the average separation of the anchors. 
In the presence of active cytoskeletal processes, the membrane 
shape fluctuations may even be increased Il23ll26l . In the mod- 
els reviewed here, the cytoskeletal elasticity is neglected since 
the relevant lateral length scales up to 50 or 100 nm are taken 
to be smaller than the average separation of the cytoskeletal 
anchors. However, the active transport of T cell receptors via 
a weak coupling to the cytoskeleton has been taken into ac- 
count in section V B and the active switching of receptors has 



been considered in section IVIII The characterization of the 
membrane elasticity by a uniform bending rigidity k is jus- 
tified on length scales larger than the molecular components 
of the membranes, i.e. on length scales larger than 5 or 10 
nm. Molecular inhomogeneities within the membranes aver- 
age out on these length scales, and the presence of anchored 
or embedded proteins leads to an increased bending rigidity, 
compared to pure lipid bilayers. 

Important length scales in the direction perpendicular to the 
membranes are the length of the receptor-ligand complexes, 
and the binding range of the receptors and ligands. The bind- 
ing range is the difference between the smallest and largest 



local membrane separation at which the molecules can bind, 
and depends on the interaction range of the molecular groups 
that stick together, the flexibility of the receptor and ligand 
molecules, and the flexibility of the membrane anchoring. In 
principle, the binding range may be measured experimentally, 
or inferred from simulations with atomistic membrane models 
in a multi-scale modeling approach. An important quantity 
is the fraction Pj, of the membranes with a local separation 
within the binding range of the receptors and ligands, see sec- 
tion [III] The membrane fraction P\, depends on the membrane 
shape fluctuations on the relevant nanoscales, and thus on the 
concentrations of the receptors and ligands, which constrain 
the shape fluctuations as bound complexes. The dependence 
of Pf, on the molecular concentrations leads to cooperative 
binding. 

As reviewed in section [ll] the integration over all possible 
positions of the receptor and ligand molecules in the parti- 
tion function of the models leads to effective adhesion poten- 
tials for the membranes. These effective adhesion potentials 
greatly simplify the characterization of the adhesion behav- 
ior In the case of long and short receptor-ligand complexes, 
for example, the effective adhesion potential allows a general 
characterization of the critical point for phase separation, see 
section [VAl If the adhesion is mediated by adsorbed particles, 
a similar integration over the degrees of freedom of these par- 
ticles leads to an effective adhesion energy that is maximal at 



an optimal particle concentration, see section VI 
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